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We investigate experimentally the early stage of the generation of waves by a tur¬ 
bulent wind at the surface of a viscous liquid. The spatio-temporal structure of 
the surface deformation is analyzed by the optical method Free Surface Synthetic 
Schlieren, which allows for time-resolved measurements with a micrometric accu¬ 
racy. Because of the high viscosity of the liquid, the flow induced by the turbulent 
wind in the liquid remains laminar, with weak surface drift velocity. Two regimes 
of deformation of the liquid-air interface are identified. In the first regime, at low 
wind speed, the surface is dominated by rapidly propagating disorganized wrin¬ 
kles, elongated in the streamwise direction, which can be interpreted as the surface 
response to the pressure fluctuations advected by the turbulent airflow. The am¬ 
plitude of these deformations increases approximately linearly with wind velocity 
and are essentially independent of the fetch (distance along the channel). Above a 
threshold in wind speed, the perturbations organize themselves spatially into quasi 
parallel waves perpendicular to the wind direction with their amplitude increasing 
downstream. In this second regime, the wave amplitude increases with wind speed 
but far more quickly than in the first regime. 

PACS numbers: 45.35.-i,47.54.-r 


I. INTRODUCTION 

Understanding the generation of surface waves under the action of wind is an old problem 
which is of primary interest for wa ve forecasting and to ev alu ate air-sea exchanges of heat, 
mass and momentum on Eartlfi^ or on natural satellites.!^ It is also important in engi¬ 
neering applications involving liquid and gas transport in pipes.!^ Despite the considerable 
literature on the subject, the physical mechanism for the onset of the first ripples at low wind 
velocity is still not fully understood. Russell,!^ as quoted by Kelvin,!^ nicely described the 
first regime where a very slight wind first destroys the perfect mirror reflection of the water 
surface, followed by a second regime for slightly larger wind where waves are observed. The 
first attempt to explain the wind-wave formation was proposed by Helmholtz and Kelvin,^!^ 
and the Kelvin-Helmholtz instability is now a paradigm for instabilities in fluid mechan¬ 
ics. However, Kelvin was aware of the discrepancy between the predicted critical wind of 
6.6 m s“^ and the commonly observed minimal wind of the order of 1 m s“^ for the first 
visible ripples on a calm sea.!^ He ascribed this discrepancy to viscous effects, which were 
not taken into account in the model. Since then, numerous attempts to better predict the 
onset of wind waves were proposed, still with limited success. 

Among the large body of literature on the subject, pioneering theoretical contributions are 
those of Phillip#^ and Miles.!^^! In an enlightening paper, PhillipJ^^ analyzed how pressure 
fluctuations in the turbulent air boundary layer could deform an otherwise inviscid fluid 
at rest. He suggested that the pressure perturbations whose size and phase velocity match 
that of the waves are selectively amplified by a resonance mechanism, and obtained a linear 
growth in time of the squared wave amplitude. The same year, MileJ^ proposed another 
mechanism based on the shear flow instability of the mean air velocity profile, ignoring 
viscosity, surface tension, drift of the liquid and turbulent fluctuations. From a temporal 
stability analysis, he showed that the boundary layer in the air is unstable if the curvature of 
the velocity profile is negative at the critical height at which air moves at the phase velocity 
of the waves, resulting in an exponential growth in time of the wave amplitude. An effort 
to classify the various instability mechanisms in parallel two-phase flow, including Miles’, is 
proposed in the review by Boomkamp and Miesen.^ 

Since then, many attempts have been made to test these predictionJi^HIZl qj- to improve 
these models,!i®^ with no definitive conclusion at the moment. While several experiments 
were devoted to determine the temporal growth of the wave after a rapid initiation of the 


2 


camera 

tr 


honej 

1 flexi 

comb 

coup 

ble 

ling 

z 

fan 





Ua^ 

V \h 


-^^ 


0 dot pattern h 


Z, = 1.5 m 

FIG. 1. Experimental setup. The wave tank and the wind tunnel are connected to the upstream 
air flow via a flexible coupling to minimize transmission of vibrations induced by the centrifugal 
fan. The surface deformations are measured by Free-Surface Synthetic Schlieren, by imaging from 
above a pattern of random dots located below the liquid tank. 


wind p^HUl other tested the amplification by wind of mechanically generated wave^^MHni or 
the wave formation by a laminar air flow.l 22 H 22 l gjjjQo boundary layers in both fluids are 
generally turbulent in the case of the air-waterJ^erfacep^HUl some authors simplified the 
problem by considering more viscous liquids 1 | 37rt 39] y^tith an airflow above a liquid more 
viscous than water, the wave onset is la rger an d, paradoxically, in better agreement with 
the inviscid Kelvin-Helmholtz prediction.EiESIini 

Rapid progresses in numerical simulations have made it possible now to address the cou¬ 
pled turbulent flows of air and water and their effect on the interface, and to access the 
pressure and stress fields hardly measurable in experiments.^^ On the experimental side, re¬ 
cent improvements in optical methods have opened the possibility to acces s exp erimentally 
the spatio-temporal structures of the waves with unprecedented resolution.EZE2l 

In the present work, we take advantage of this technical improvement to analyze the 
early stage of wave formation at the surface of a viscous liquid. Surface deformations are 
measured with a vertical resolution better than one micrometer using Free-surface Synthetic 
Schlieren,!^ a time-resolved optical method based on the refraction of a pattern located below 
the fluid interface. Working with a viscous liquid has two advantages: first, the flow in the 
liquid remains laminar and essentially unidirectional with a limited surface drift; second, 
the perturbations of the interface that are not amplified by an instability mechanism are 
rapidly damped, so the surface deformations at low wind velocity are expected to be the 
local response in space and time to the instantaneous pressure fluctuations in the air. Our 
results clearly exhibit two wave regimes: (i) at low wind velocity, small disordered surface 
deformations that we call ’’wrinkles” first appear, elongated in the streamwise direction, 
with amplitude growing slowly with the wind velocity but with no significant evolution with 
fetch (the distance upon which the air blows on the liquid); (ii) above a well defined wind 
velocity, a regular pattern of gravity-capillary waves appears, with crests normal to the wind 
direction and amplitude rapidly increasing with wind velocity and fetch. 


II. EXPERIMENTAL SET-UP 
A. Liquid tank and wind tunnel 

The experimental set-up is sketched in Fig. It is composed of a fully transparent 
Plexiglas rectangular tank of length L = 1.5 m, width W = 296 mm, and depth h = 35 mm, 
fitted to the bottom of a horizontal channel of rectangular cross-section. The channel width 
is identical to that of the tank, its height is H = 105 mm, with two horizontal floors of 
length 26 cm before and after the tank. The tank is filled with a water-glycerol mixture, 
such that the surface of the liquid precisely coincides with the bottom of the wind tunnel. 

Air is injected upstream by a centrifugal fan through a honeycomb and a convergent (ratio 
2.4 in the vertical direction). To minimize transmission of vibrations induced by the fan, the 
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wind-tunnel is mounted on a heavy granite table and connected to the upstream channel via 
a flexible coupling. Residual vibrations induce surface deformations less than 1 /rm. The 
wind velocity Ua, measured at the center of the outlet of the wind tunnel with a hot-wire 
anemometer, can be adjusted in the range f — 10 m s“^. We define x in the streamwise 
direction (fetch), y in the spanwise direction and z in the vertical direction. The origin 
( 0 , 0 , 0 ) is located at the free surface at fetch 0 , at mid-distance between the lateral walls. 

The tank is filled with a mixture of 80% glycerol and 20% water, of density p = 1.20 x 
10^ kg m“^ at 25°C (the room temperature being regulated to this temperature). Kinematic 
viscosity, measured with a low shear rheometer, is v = rj/p = 30 x 10“® m^ s“^ at this 
temperature. The water-glycerol mixture is extremely sensitive to surface contamination, 
which may induce strong surface tension gradients and alter both the mean flow in the liquid 
and the generation of waves.^ To overcome this problem, we let the wind blow for a few 
minutes, and we remove the contaminated part of the surface liquid by collecting it at the 
end of the tank. The procedure is repeated frequently, and in normal operating conditions 
the surface of the liquid remains clean over most of the liquid bath, with less than 30 cm 
of polluted surface remaining at the end of the tank. Surface tension of the clean mixture, 
measured with a Wilhelmy plate tensiometer, is 7 = 60 ± 5 mN m“^, and the capillary 
wavelength is Ac = 2Try^j/pg ~ 14.2 mm. The dispersion relation for free surface waves 
propagating in an inviscid liquid at rest is 

= ^gk + tanh{kh), ( 1 ) 

where w and k are the angular frequency and wave number. Finite depth effects can be ne¬ 
glected in the present experiment: the depth correction factor, tanh(fch), is larger than 0.98 
for wavelength smaller than 90 mm. In spite of the large viscosity used in our exp erime nts, 
viscous correction to the inviscid dispersion relation can be also neglected hereP^^^ the 
phase velocity matches the inviscid prediction to better than 10 “^ for the waves observed at 
onset (A ~ 30 mm). On the other hand, this large viscosity induces a strong attenuation of 
the waves. For the wave tank geometry and the typical wavelengths considered here, friction 
with the bottom and side walls is negligihle, and the attenuation length for free waves is 
governed by the dissipation in the bulk,HSl = Cg/(2j^fc^), with Cg{k) the group velocity. 
For A ~ 30 mm the attenuation length is Ly ~ 60 mm, indicating that a free disturbance at 
this wavelength cannot propagate over a distance much larger than a few wavelengths. As 
a consequence, although the tank is of limited size, reflections on the walls or at the end of 
the tank can be neglected in our experiment. 


B. Wind profile 


The velocity profile in the air U{z), measured using hot-wire anemometry, is shown in 
Fig. i for a wind velocity Ua = 3.9 m s ^ at fetch x = 20, 500 et 1000 mm. The hot-wire 
(Dantec Dynamics 55P01) is 5 /rm in diameter with an active length of 1.25 mm, and is 
mounted on a sliding arm to allow vertical motion with a 0.1 mm accuracy. The velocity 
profiles show the development of the boundary layer along the channel: the thickness do, 99 , 
defined as the distance from the surface at which the mean velocity is 0.9917a, increases 
nearly linearly, from 12.6 mm at a; = 0 to 32 mm at a; = 1.0 m (slope of order of 2%). 
The fact that do. 99 (a;) approaches the channel half-height i7/2 ~ 52 mm at the end of the 
channel indicates that the flow becomes fully developed there. 

The evolution of the friction velocity u*{x) along the channel can be ob tained by fitting 
the velocity profiles for z < do. 99 (a;) with the classical logarithmic law, h^ l ^^ l 


U{x, z) 
u*{x) 



+ C, 


( 2 ) 


with K ~ 0.4 the Karman constant, C = 5, and Sy(x) = i'a/u*{x) the thickness of the viscous 
sublayer. We find u* to slightly decrease with fetch: for Ua = 3.9 m s“^, u* decreases from 
0.22 m s“^ at a: ~ 0 down to 0.17 m s“^ at a: = 1 m. Accordingly, dt,(a:) slightly increases 
with fetch, from 0.07 to 0.09 mm. 

The procedure is repeated for different wind velocities at a fixed fetch, xq = 500 mm. 
Measurements are restricted to C/q < 6 m s“^, when the surface deformations remain weak 
(less than 10 ^m), because the hot-wire could not be positioned too close to the liquid. We 
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FIG. 2. Mean velocity profiles U ( 2 ) in the air for Ua = 3.9 m at fetch x = 20, 500 and 1000 mm 
and y = 0. The dotted line shows the streamwise development of the 99% boundary-layer thickness, 
So. 99 ix) — 12.6 mm + 0.02a;. The curves for 2 < 5Q.gg(x) show the fit with the logarithmic law 11- 

find that in this range u* is almost proportional to Ua , u*{xo) — 0.0517a (see inset in Fig.jl. 
The corresponding half-height channel Reynolds number at this fetch, Rcr = Hu*/2va, 
varies in the range 160 — 1000, and the thickness of the viscous sublayer 5y decreases from 
0.3 to 0.05 mm when Ua increases from 1 to 6 m s“^. Since the flow in the viscous sublayer is 
essentially laminar up to z ~ lO^i,, which is comfortably larger than any surface deformation 
over this range of velocity, we can consider the air flow to be close to a canonical turbulent 
boundary over a no-slip flat wall, at least for a wind velocity up to 6 m s“^. This does not 
hold for larger wind veloc ity, for which the roughness induced by the waves decreases the 
value of C in Eq. 


C. Flow in the liquid tank 


The shear stress induced by the wind at the interface drives a drift flow in the liquid. 
Since the tank is closed, this drift is compensated by a back-flow at the bottom of the tank, 
and a stationary state is reached after a few minutes. We have measured the mean velocity 
profile in the tank using Particule Image Velocity (PIV) in vertical planes {x,z) (Fig. [^. 
Except at small fetch (on a distance of the order of the liquid height) and over the last 30 cm 
of the tank (where surface contamination cannot be avoided), the velocity profiles are found 
nearly homogeneous in x and y. The velocity profiles are well described by the parabolic 
law, solution of the stationary Stokes problem 

u{x,z) = Us{x) (1 ++ , (3) 

for —h < z < 0, where Us{x) = u{x,z = 0) is the surface velocity. For 17o = 4 m s“^the 
surface velocity is of order of 1 cm s“^, which leads to a Reynolds number Re = Ugh/v ~ 10 
in the liquid. This drift velocity is in agreement with measurements at small Reynolds 
number,E3 but is muc h smaller than the 2-3% of wind velocity typically found in classical 
air-water experiments.^^^EHIZlSnixhe small surface velocity here is expected to have negligible 
effect on the dispersion relation Q: Lilly (see appendix of Hidy and Plat^^ shows that 
the correction to the phase velocity for this parabolic profile is 2Ua/kh, which is only 2 % of 
the phase velocity for the most unstable wavelength (A ~ 30 mm). 

Because of the development of the boundary layer and the resulting decreasing friction 
velocity, the surface velocity U^ decreases slightly along the tank. For Ua = 4 m s“^, Ua{x) 
decreases from 1.5 to 1.1 cm s“^. Measuring Us{x) provides another way to determine u*{x): 
using the continuity of the stress at the interface, one has a{x) = paU*'^{x) = r]du/dz{z = 0 ), 
yielding 



u 


( 4 ) 
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FIG. 3. Mean velocity profiles in the liquid measured by PIV at fetch x = 400 mm for various wind 
velocities Ua- The profiles are averaged in time and in the streamwise direction over Ax = 100 mm. 
The continuous line for the largest value of Ua shows the quadratic profile 0- Inset: friction 
velocity u*, deduced from the mean profile in the airflow at xo = 500 mm (squares) and deduced 
from the shear stress at the liquid surface (Eq. (circles), as a function of the wind velocity Ua- 
The continuous line is a fit by u* = 0.05C/a. 


The friction velocity u* measured by PIV in the liquid with this method is in excellent 
agreement with the one measured in air with the hot-wire at Xq = 500 mm (inset of Fig.j^. 
For simplicity, the ratio of u*/Ua is taken in the following as constant and equal to 0.05 for 
all fetches. By comp arison, the ratio of u*/Ua is generally found of order of 3% in air-water 
experimentsj ^^ l ^^ ^^ with weak dependence on the wind velocity.!^ 

Note that the shear stress at the liquid surface and at the lateral and upper walls must 
be balanced by a small longitudinal pressure gradient Ap/L in the air along the channel. 
This pressure gradient introduces a complication in the setup: the liquid surface becomes 
slightly tilted, with the inlet liquid height below the outlet height (this is analogous to the 
’wind tide’ effect observed on lake^^. Assuming equal stress a on the liquid surface and on 
the solid walls, this pressure gradient writes Ap/L ~ 2cr(l/VF -I- 1/H), with W and H the 
channel width and height. For a wind velocity C/a = 4 m s“^, the pressure drop along the 
tank is Ap ~ 2 Pa, which results in a hydrostatic height difference between the two ends 
of the tank of Ap/pg ~ 0.2 mm, in good agreement with our measurement. We observed 
that the resulting backward facing step that could appear at a: = 0 increases the turbulent 
fluctuations and significantly enhances the wave amplitude at small fetch by typically a 
factor of 2. It is therefore critical to maintain the liquid level at x = 0 by carefully tilting 
the channel. We achieve a leveling of the liquid at x = 0 better than 20 /im by using the 
tangential reflexion of a laser sheet intersecting the upstream plate and the liquid surface. 


D. Surface deformation measurement 

We measure the surface deformation of the liquid using the Free Surface Synthetic 
Schlieren (FS-SS) method.l^ This optical method is based on the analysis of the refracted 
image of a pattern visualized through the interface. A random dot pattern located below 
the liquid tank is imaged by a fast camera located above the channel, with a field of view of 
390 X 280 mm. A reference image is taken when the liquid surface is flat (zero wind), and 
the apparent displacement field 5r between this reference image and the distorted image in 
the presence of waves is computed using an image correlation algorithm. Integration of this 
displacement held gives the height held C{x,y,t) (see examples in Fig.|^. 

Measurements are performed at three fetches, corresponding to the the flrst three quarters 
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FIG. 4. Instantaneous surface height C,{x, y) measured by FS-SS centered at intermediate fetch x = 
570 mm, at increasing wind velocities, (a) Ua = 3.2 m s“^, showing small-amplitude disorganized 
wrinkles elongated in the streamwise direction (^rms = 0.0032 mm), (b) Ua = 5.9 m s“^, showing 
a combination of streamwise wrinkles and spanwise waves ((^rms = 0.009 mm), (c) Ua = 7.0 m s”'^, 
showing well-defined spanwise waves of mean wavelength A = 35 mm (Crms = 0.12 mm), (d) 
Ua = 7.8 m s“^, showing large-amplitude waves of mean wavelength A = 44 mm with increasing 
disorder (^rma = 0.6 mm). Note the change of scale in the color map. 


of the tank with a small overlap: x € [10,400] mm; x G [370, 760] mm; x G [700,1090] mm. 
No measurements are performed in the last quarter of the channel because of possible surface 
contamination. The distance between the random dot pattern and the liquid surface sets 
the sensitivity of the measurement, and is chosen according to the typical wave amplitude. 
We chose a distance of 29 cm for waves of weak amplitude (of order of 1 — 10 f^m), and 
6 cm for waves of large amplitude (up to 1 mm). For wave amplitude larger than a few 
millimeters, the FF-SS method no longer applies: crossing of light rays appear below waves 
of large curvature (caustics), which prevents the measurement of the apparent displacement 
field. The horizontal resolution is 3 mm, and the vertical resolution of order of 1% of 
the wave amplitude. Acquisitions of 2 s at 200 Hz are performed for time-resolved wave 
reconstruction, and 100 s at 10 Hz to ensure good statistical convergence of the root mean 
square of the wave amplitude. 


III. RESULTS 

A. Amplitude versus velocity 

Figure shows four snapshots of the surface deformation at increasing wind velocity, 
between Ua = 3.2 and 7.8 m s“^, at intermediate fetch x G [370,760] mm. At small Ua 
the wave pattern shows rapidly moving disorganized wrinkles of weak amplitude, of order 
10 /xm, elongated in the streamwise direction (Fig. |^). As the wind velocity is increased, 
noisy spanwise crests, normal to the wind direction, gradually appear in addition to the 
streamwise wrinkles (Fig. &)• The amplitude of these spanwise crests rapidly increases 
with velocity and becomes much larger than the amplitude of the streamwise wrinkles for 
wind velocity in the range 6 — 7 m s“^. At Ua = 7.0 m s“^(Fig. |^) the surface field 
is dominated by a regular wave pattern of typical amplitude 0.2 mm, with a well defined 
wavelength in the streamwise direction. The wave crests are not strictly normal to the wind. 
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FIG. 5. Root mean square of the surface height (^rms as a function of wind velocity Ua- The data 
are averaged over the measurement windows centered at two values of the fetch x. The vertical 
arrows show the velocities corresponding to the four snapshots in Fig.|^ The continuous line shows 
the linear fit = cdJa, with a = 10~® s. 


but rather show a dislocation near the center line. This may be due to a slight convergence 
of the turbulent air flow close to the free surface towards the walls, which is unavoidable 
for a turbulent channel flow in a rectangular geometry (secondary flow of Prandtl’s second 
kincP^. As the wind speed is further increased, the disorder of the wave pattern increases 
(Fig. gi), with more dislocations and larger typical wavelength and amplitude. All these 
patterns are quite similar to those reported by Lin et a/.^ from direct numerical simulation 
of temporally growing waves with periodic boundary conditions. 

The evolution from the disorganized longitudinal wrinkles to the well-defined transverse 
waves as the wind velocity is increased is evident from the root mean square of the defor¬ 
mation amplitude, 

Crms = 

where the brackets are both temporal average and spatial average over the field of view. 
This quantity, plotted as a function of the wind velocity in Fig. for two values of the 
fetch X, clearly exhibits the two regimes: at small air velocity, when the surface deformation 
is dominated by the longitudinal wrinkles, the wave height slowly increases with the wind 
velocity, but beyond a threshold of order of 6 m s“^ the increase becomes much sharper and 
fetch dependent: the wave amplitude grows by a factor of 100 for Ua increasing between 6 and 
8 m s“^. A similar transition in the wave amplitude is also reported in air-water experiments 
by Kahma and DonelarP^ and Caulliez et al.^^ At the largest velocity. C/a — 8 m s“^, the 
sharp increase of the wave amplitude apparently starts to saturate. This wind velocity 
represents an upper limit for the FS-SS measurements because of the caustics induced by 
the strong wave curvature. 

In the wrinkle regime (C/q < 6 m s“^), the wave height is almost independent of the fetch, 
and is approximately proportional to the wind velocity: Crms — ctUa, with a = 10“® s. 
This independence of x suggests that the wrinkles can be simply viewed as an imprint 
on the free surface of the turbulent fluctuations in the airflow. Relating quantitatively 
the height fluctuations to the pressure fluctuations is however a difficult task. A simple 
estimate, assuming an instantaneous hydrostatic response of the liquid interface (i.e., ne¬ 
glecting viscous and capillary effects) would yield ^rms — Prms/pg- The pressure fluctu ation 
at the wall in a fully developed turbulent channel is well described by the empirical lavEsMl 
Prms = /(-Rer)PaU*^, with f {Rcr) = (2.60 ln(i?eT-) —11.25). In the range Ua — 1 — 6 m s“^, 
one has Rcr — 160 — 1000 , which (neglecting the logarithmic variation over this range and 











FIG. 6. Wave amplitude ^rms (averaged over the spanwise coordinate y and time) as a function of 
the fetch x for various wind speeds Ua- In the wrinkle regime {Ua < 6.3 m s“^) Crins(®) is almost 
constant, whereas it increases with fetch in the wave regime. 


taking u* ~ 0.05C/a) yields Prms — O.OOdpaU^, and hence Crms — 0.3 — 20 /rm. Although 
the order of magnitude is consistent with Fig. the predicted scaling (Crms oc U^) is not 
compatible with the data, suggesting that the viscous time response of the liquid must be 
accounted for to describe the observed trend Crms oc Ua- 


B. Spatial growth rate 

Contrarily to the wrinkles, which are almost independent of the fetch x, the amplitude of 
the transverse waves strongly increases with x, as shown in Fig.j^ The rms amplitude Crms(a;) 
is computed here using an average over y and time only. The spatial growth is approximately 
exponential at small fetch (x < 400 mm), as expected for a convective supercritical instability 
in an open flow.l^ At larger fetch nonlinear effects come into play, resulting in a weaker 
growth of the wave amplitude. 

The spatial growth rate /3 can be estimated in the initial exponential growth regime 
(a; < 400 mm) by fitting the squared amplitude as CSns(^) exp(/3a:). The growth rate ,0, 
plotted in Fig. as a function of the wind velocity, allows to accurately define the onset of 
the wave growth: one has /3 ~ 0 for Ua < 6.3 m s“^, and a linear increase at larger Ua, 
which can be fitted by 


P ~ b{Ua - Ua) 


( 5 ) 


with Ua — 6.3±0.1 m s“^ and b ~ 11.6±0.8 s m~^. Interestingly, in the small range of fetc h 
where /3 is computed, the waves are nearly monochromatic, with A ~ 30 mm (see Sec. Ill C). 
This indicates that the growth rate measured here, although computed from the total wave 
amplitude, corresponds essentially to the growth rate of the most unstable wavelength. 

We note that the velocity threshold Ua — 6.3 m s“^ turns out to be close to the (inviscid) 
Kelvin-Helmholtz prediction. This agreement, first noted by Franci^^for a viscous fluid, is 
however coincidental since the threshold depends on viscosity.l^ 

An interesting question is whether the wrinkles at low wind velocity can be considered as 
the seed noise for the exponential growth of the waves at larger velocity. Figure [^indicates 
that this is apparently not the case: for Ua > 6.3 m s“^, the initial wave amplitude (^rms 
extrapolated at a: = 0 increases with Ua much more rapidly than the amplitude of the 
wrinkles; Crms(a: = 0) grows from 8 to 30/im for Ua increasing from 6.3 to 7.7 m s“^ only. 
This suggests that the wrinkles are not necessary for the growth of the waves. Instead, the 
seed noise for the waves probably results from the surface disturbance induced by the sudden 
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FIG. 7. Spatial growth rate jS, measured at small fetch {x < 400 mm), as a function of wind speed. 
The continuous line corresponds to the linear fit ill- 


change in the boundary condition from no-slip to free-slip at a; = 0. Accordingly, the rms 
amplitude can be described as the sum of the wrinkle amplitude (linearly increasing with 
Ua) and the wave amplitude (exponentially increasing with Ua), 

Crms ^ alia + Cn{Ua) BXp [6 {Ua - Uc) x] , (6) 

with Cn{Ua) the amplitude of the noise at zero fetch. 

To provide comparison with other experiments and theoretical results, it is interesting to 
express Eq. (§ in terms of a temporal growth rate. In the frame moving with the group 
velocity, this temporal growth rate writes (it = Cg/3, with Cg ~ 0.16 m s“^ for the most 
unstable wavelength A ~ 30 mm. In terms of the friction velocity u* ~ O.OSt/a, Eq. ([^ 
writes 


Pt ^ (36±8)(m* - 0.31) 

(in s“^). Not su rprisingly, these values are smaller than the ones reported in air-water 
experimentsby a factor of order of 3, suggesting that the wave growth is weakened 
by the viscosity of the liquid. 


C. Spatial structures 


To characterize the spatial structure of the wrinkles and waves, we introduce the two-point 
correlation 


N _ (C(x,t)C(x-fr,t)) 

where (•) is a spatial and temporal average. The correlation in the streamwise direction 
(r = Tx^x) is plotted in Eig. [^for the four wind velocities corresponding to the snapshots 
in Eig. 1^ The monotonic decay of C{rx) at small wind velocity is a signature of the 
disordered deformation pattern in the wrinkle regime, whereas the oscillations at larger 
velocity characterize the onset of waves. Interestingly, these oscillations are clearly visible 
even at Ua = 5.9 m s“^, confirming that the transverse waves are already present in the 
deformation field significantly before the critical velocity [4 — 6.3 m s“^(see Eig. |^). This 
indicates that the transition between the wrinkles and the waves is not sharp: both structures 
can be found with different relative amplitude over a significant range of wind velocity. 
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FIG. 8. Spatial correlation function in the streamwise direction, C{rx), at four wind velocities for 
X G [370, 760] mm. 


The smooth transition between wrinkles and waves can be further characterized by com¬ 
puting the correlation length in the direction [i = x,y), which we define as 6 times 
the first value of satisfying C{ri) = 1/2. This definition is chosen so that coincides 
with the wavelength for a monochromatic wave propagating in the direction e^. Although 
no wavelength can be defined for the disorganized wrinkles, A^, and Ay provide estimates 
for the characteristic distance between wrinkles in the streamwise and spanwise directions. 

The correlation lengths A^, and Ay are shown in Fig.|^as functions of the wind velocity. At 
very low wind (~ 1 m s“^), both lengths are of the same order, A^, ~ Ay ~ 250 mm. Between 
1.5 and 5 m s“^, the surface deformations are mostly in the streamwise direction [A^/Ay ~ 
3), whereas at larger velocity they are essentially in the spanwise direction {A^/Ay ~ 0.15). 
Note that pure monochromatic waves in the x direction would yield C{ry) = 1 and hence 
Ay = oo; the saturation of Ay close to the channel width {W = 296 mm) observed at large Ua 
is a signature of the dislocation existing near the center line y = 0 and visible in Fig. ic,d). 

It is worth noting that the increase of Ay and the decrease of A^, start at a wind velocity 
L/j ~ 5 m s“^ which is significantly lower than the critical velocity L/, ~ 6.3 m s“^, 
confirming that transverse waves are present well before their amplification threshold. This 
coexistence of wrinkles and waves at Ua < Uc suggests the following picture: below the onset, 
waves are locally excited by the wrinkles, but they are exponentially damped (/3 < 0). The 
surface field can therefore be described as the sum of a large number of spatially decaying 
transverse waves, locally excited by the randomly distributed wrinkles generated by the 
pressure fluctuations in the boundary layer. Since the amplitude of the wrinkles is essentially 
independent of x, the resulting mixture of wrinkles and decaying transverse waves is also 
independent of x, leading to the apparent growth rate /3 = 0 of Fig. In other words, 
the expected negative growth rate below the onset is hidden by the spatial average over 
randomly distributed decaying waves, and cannot be inferred from the observed constant 
deformation amplitude for Ua < Uc- 

If the elongated wrinkles at low velocity are traces of the pressure fluctuations in the 
boundary layers, we expect a relationship between their characteristic dimensions. The 
geometrical and statistical properties of the pressure fluctuations in a turbulent channel 
cannot be obtained experimentally, but are available from numerical simulations. We refer 
here to data from Jimenez and Hoya^^ at Rcr up to 2000. The intensity of the pressure 
fluctuations increases logarithmically from the center of the channel down to z ~ 30i5„, 
and then remains essentially constant in the viscous sublayer down to z = 0 (see their 
Fig. 8b). In the thin region where the pressure fluctuations are maximum, 0 < z < 30i5„, 
the characteristic dimensions of the pressure structures in the {x,y) planes normal to the 
wall are nearly equal, — £y — 160Jt,, and are hence decreasing with increasing wind 
velocity. These features are indeed compatible with the correlation lengths in Fig. at 
very low velocity the wrinkles are nearly isotropic, A^, ~ Ay ~ 250 mm. As Ua increases up 
to 3 m s“^, the spanwise correlation length decreases as Ay ~ 550i5^, similarly to the width 
iy of the pressure fluctuations. This decrease is however not observed for the streamwise 
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FIG. 9. Streamwise and spanwise correlation lengths A^; (□) and Ay (a), averaged over x € 
[700,1090] mm, as a function of the wind velocity Ua- The vertical dashed line at Ua = 6.3 m s“^ in¬ 
dicates the onset of the wave growth. The continuous line is a fit Ay ~ 5505^, with Sv = Vaju* the 
thickness of the viscous sublayer. 


correlation length A^,, which may result from the viscous time response of the surface behind 
a moving pressure perturbation. 


D. Spatiotemporal dynamics 


In order to confirm the relation between the surface wrinkles and the pressure fluctuations 
traveling in the boundary layer, we now turn to a spatio-temporal description of the surface 
deformation. We show in Fig. 10 spatio-temporal diagrams (left) and two-point two-time 
correlation (right) at increasing wind velocity. The spatio-temporal diagrams are constructed 
by plotting the surface deformation ({x,y,t) in the plane {x,t) along the center line y = 0- 
The oblique lines in these diagrams indicate the characteristic velocity of the deformation 
patterns. The spatio-temporal correlation in the streamwise direction is dehned as 


Cirx,T) = 


(C(x + ry;ey,,t-f r)C(x,t)) 

(C(x,t)2) 


( 7 ) 


where (•) is a spatial and temporal average. For statistically stationary and homogeneous 
deformations, one has C{—rx, —t) = C(rx,T), so only the positive time domain is shown. 

At small wind velocity {Ua = 3.2 m s“^) the surface deformation shows rapidly propagat¬ 
ing disorganized structures, with life time of order of their transit time [Fig. [Ioj]a)]. Their 
characteristic velocity is distributed over a large range, resulting in a broad correlation in 
Fig.[^b). 

In the mixed wrinkle-wave regime {Ua = 5.9 m s“^), slow wave packets with well defined 
velocity appear, embedded in a sea of rapid disorganized fluctuations [Fig. [To|[c)]. These 
slow wave packets confirm the picture of transverse waves locally excited by the wrinkles 
but rapidly damped because of their negative growth rate. The wavelength and the phase 
velocity of these evanescent waves can be inferred from the corresponding spatio-temporal 
correlation [Fig. [To|(d)]. 

At larger wind velocity {Ua — 7.0 m s“^), the surface deformation becomes dominated 
by spatially growing transverse waves (/3 > 0), resulting in well dehned oblique lines in the 
spatio-temporal diagram [Fig. |10[ e)] and a marked spatial and temporal periodicity in the 
correlation [Fig. |10| [f)]. These transverse waves, however, are never strictly monochromatic: 
wave packets are still visible, delimited by boundaries propagating at the group velocity. 
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FIG. 10. (a,c,e,g) Spatio-temporal diagrams ^{x,t) taken along the line y = 0. Same velocities and 
scales as in Fig. (b,d,f,h) Longitudinal spatio-temporal correlation C{rx,T) [Eq. 01-, Colormap 
is [—1,1] from blue to red. The black dashed line shows the convection velocity T4onv (|^. In (f), 
the two white dashed lines show the group velocity Cg = 0.16 m s“^ corresponding to the observed 
wavelength X — 37 mm, which delimit the wave packets. 


For this wind velocity and fetch, the local wavelength is A = 37 mm, for which the predicted 
phase and group velocities are c = 0.26 m s“^ and Cg = 0.16 m s“^, in good agreement 
with the observed slopes (black and white dashed lines, respectively) in Fig. ©f)- 

Finally, at even larger wind velocity {Ua = 7.9 ms the phase velocity and wavelength 
increase slightly, and becomes again broadly distributed. Accordingly, the spatial and tem¬ 
poral periodicity weakens in Fig. |10[ h). 

For monochromatic waves propagating in the x direction, one has C(rx, t) = cos{krx — ojt), 
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FIG. 11. Convection velocity Konv = A^j/r, averaged over x G [700,1090] mm, as a function of the 
wind velocity Ua- The vertical dashed line at Uc = 6.3 m s“^ indicates the onset of wave growth, 
13 = 0. The continuous line shows Korr = 0.61/a. The dashed horizontal line shows the local phase 
velocity c = 0.27 m s“^ corresponding to the local wavelength A = 44 mm. 


SO the correlation is 1 along characteristic lines parallel to r^/t = c, where c = oj/k is the 
phase velocity. For a (non-monochromatic) propagating pattern, the correlation is weaker 
but remains maximum along a line Vx/t given by the characteristic velocity of the pattern. 
We therefore define the convection velocity as 


K 

V — — 

^conv — 5 

r 


( 8 ) 


where A^, and r are defined as the first value such that C{h.x/0^ 0 ) = 1/2 and C( 0 , t/6) = 1 / 2 . 
Here again, the factor 6 is chosen such that A^, and r correspond to the wavelength and 
period for monochromatic waves; in this case, V/onv is simply the phase velocity. This 
convection velocity is shown as black dashed lines in the spatio-temporal correlation in 
Fig. [T^ and is plotted in Fig. El as a function of the wind velocity Ua- At very small 
velocity {Ua < 2 m s“^), V/onv first increases, following approximately V/onv — O-dUa- This 
value is remarkably similar to the convection velocity of the pressure fluctuations found in 
turbulent boundary layers.^ This is consistent with the fact that, in this small range of Ua, 
the decrease of the characteristic width of the wrinkles Ay follows the expected decrease 
of the size of the pressure fluctuations in the boundary layer (Fig. |9|. Such convection 
velocity 14onv — 0.617a can be recovered as follows: in a turbulent boundary layer the 
pressure fluctuations are maximum at ZmlSy — 20 — 50 (see, e.g., Jimenez and Hoyas^ and, 
according to the logarithmic law the mean velocity at this height is u{zm) — (13± l)u*. 
Using u* ~ 0.0517a in the present experiment, this yields u{zm) — (0.65 ± 0.05)17a. This 
suggests that the surface response to the traveling pressure fluctuations is essentially local 
and instantaneous up to 17a — 2 m s“^. 

For t7a > 2 m s“^, the convection velocity departs from the linear growth 0.617a, and 
saturates at Wonv — 1.2 m s“^. This saturation probably results from the viscous damping, 
which prevents an instantaneous response of the surface deformation to too rapidly prop¬ 
agating pressure disturbances. As the wind velocity is further increased, the convection 
velocity decreases down to c ~ 0.27 m s“^, which coincides with the expected phase velocity 
of free waves of the observed wavelength. This gradual decrease does not correspond to a 
slowing of the wrinkles, but rather to an average over the rapid wrinkles propagating at 
velocities of the order of 1.2 m s“^, which dominate the surface at low Ua, and the slow 
transverse waves at velocity of the order of 0.27 m s“^, which dominate the surface at high 
Ua- Finally, for t7a > 7 m s“^, the convection velocity increases again, which is consistent 
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with the increase of the wavelength in Fig. These increases do not necessarily occur for 
shorter fetches, implying that the wave properties change with fetch at high wind veloc¬ 
ity. Such non linear effect presents strong similarities with the wavenumber and frequency 
downshift observed for wind-waves on sea, and will be investigated in future work.” 


IV. CONCLUSION 

In this paper, we explored the spatio-temporal properties of the first surface deforma¬ 
tions induced by a turbulent wind on a viscous fluid. New insight into the wave genera¬ 
tion mechanism is gained from spatio-temporal correlations computed from high resolution 
time-resolved measurements of the surface deformation field. At low wind velocity, rapidly 
propagating disordered wrinkles of very small amplitude are observed, resulting from the 
response of the surface to the traveling pressure fluctuations in the turbulent boundary layer. 
Above a critical wind velocity Uc, we observe the growth of well defined propagating waves, 
with growth rates compatible with a convective supercritical instability. Interestingly, an 
intermediate regime with spatially damped waves locally excited by the wrinkles is observed 
below Uc, resulting in a smooth evolution of the characteristic lengths and velocity as the 
wind speed is increased. Above the onset t/c, the seed noise for the growth of the waves 
is apparently not governed by the wrinkles, but rather by the perturbations at the inlet 
boundary condition at zero fetch. 

Using a liquid of large viscosity yields considerable simplification of the general problem 
of wave generation by wind. Although some of the present results may be relevant to the 
more complex air-water configuration, other are certainly specific to the large viscosity of 
the liquid. In particular, the wrinkles observed at very low wind velocity (Ua < 2 m s“^) are 
compatible with a local and instantaneous response of the surface to the pressure fluctuations 
traveling in the boundary layer. This simple property is not expected to hold for liquids 
of lower viscosity such as water, for which the surface deformation at a given point results 
from the superposition of the disturbances emitted previously from all the surface. New 
experiments with varying viscosity are necessary to gain better insight into the intricate 
relation between the turbulent pressure field and the surface response below the onset of the 
wave growth, and to characterize the spatial evolution of the waves above onset. 


ACKNOWLEDGMENTS 

We are grateful to H. Branger, C. Clanet, P. Clark di Leoni, B. Gallet, J. Jimenez, 
O Naraigh, and P. Spelt for fruitful discussions. We acknowledge A. Aubertin, L. Auffray, 
C. Borget, and R. Pidoux for the design and set-up of the experiment. This work is supported 
by RTRA ’’Triangle de la Physique”. F.M. thanks the Institut Universitaire de France. 

^P. H. LeBlond and L. A. Mysak. Waves in the Ocean. Elsevier, 1981. 

^P. Janssen. The interaction of ocean waves and wind. Cambridge University Press, 2004. 

®AG Hayes, RD Lorenz, MA Donelan, M Manga, JI Lunine, T Schneider, MP Lamb, JM Mitchell, WW Fis¬ 
cher, SD Graves, et al. Wind driven capillary-gravity waves on Titan’s lakes: Hard to detect or non¬ 
existent? Icarus, 225:403—412, 2013. 

■^J. W. Barnes, G. Sotin, J. M. Soderblom, R. H. Brown, A. G. Hayes, M. Donelan, S. Rodriguez, 
S. Le Mouelic, K. H. Baines, and T. B. McCord. Cassini/vims observes rough surfaces on titan’s punga 
mare in specular reflection. Planetary Science, 3(1):1-17, 2014. 

®G. Hewitt. Annular two-phase flow. Elsevier, 2013. 

®J Scott Russell. On waves. In Report of fourteenth meeting of the British Association for the Advancement 
of Science, York, pages 311—390, 1844. 

^W. Thomson. On the waves produced by a single impulse in water of any depth, or in a dispersive medium. 

Proceedings of the Royal Society of London, 42:80—83, 1887. 

®H. L. von Helmholtz. On discontinuous movements of fluids. Philos. Mag., 36:337, 1868. 

®0. Darrigol. Worlds of Flow: A History of Hydrodynamics from the Bemoullis to Prandtl. Oxford 
University, 2005. 

^'^O. M. Phillips. On the generation of waves by turbulent wind. J. Fluid Mech., 2(05):417-445, 1957. 

J. W. Miles. On the generation of surface waves by shear flows. J. Fluid Mech., 3:185-204, 1957. 

^^P. Boomkamp and R. Miesen. Classification of instabilities in parallel two-phase flow. International 
Journal of Multiphase Flow, 22:67—88, 1996. 

^®E. J. Plate, P. C. Chang, and G. M. Hidy. Experiments on the generation of small water waves by wind. 
J. Fluid Mech., 35(4):625-656, 1969. 

^■^K. Kahma and M. A. Donelan. A laboratory study of the minimum wind speed for wind wave generation. 
J. Fluid Mech., 192:339-364, 1988. 


15 


Liberzon and L. Shemer. Experimental study of the initial stages of wind waves’ spatial evolution. J. 
Fluid Mech., 681:462-498, 2011. 

Grare, W. Peirson, H. Branger, J. Walker, J-P. Giovanangeli, and V. Makin. Growth and dissipation 
of wind-forced, deep-water waves. J. Fluid Mech., 722:5—50, 2013. 

S. Hristov, S. D. Miller, and C. A. Friehe. Dynamical coupling of wind and ocean waves through 
wave-induced air flow. Nature, 422(6927):55—58, 2003. 

Katsis and T. R. Akylas. Wind-generated surface waves on a viscous fluid. Journal of applied mechanics, 
52(1):208-212, 1985. 

Teixeira and S.E. Belcher. On the initiation of surface waves by turbulent shear flow. Dynamics of 
atmospheres and oceans, 41(l):l-27, 2006. 

^^G. R. Valenzuela. The growth of gravity-capillary waves in a coupled shear flow. J. Fluid Mech., 
76(02) :229-250, 1976. 

van Gastel, P. Janssen, and G. J. Komen. On phase velocity and growth rate of wind-induced gravity¬ 
capillary waves. J. Fluid Mech., 161:199-216, 1985. 

^^W. R. Young and G. L. Wolfe. Generation of surface waves by shear-flow instability. J. Fluid Mech., 
739:276-307, 2014. 

Mitsuyasu and K. Rikiishi. The growth of duration-limited wind waves. J. Fluid Mech., 85(04) :705-730, 
1978. 

Kawai. Generation of initial wavelets by instability of a coupled shear flow and their evolution to wind 
waves. J. Fluid Mech., 93(4):661-703, 1979. 

Veron and W. K. Melville. Experiments on the stability and transition of wind-driven water surfaces. 
J. Fluid Mech., 446(10):25-65, 2001. 

B. Bole and E. Y. Hsu. Response of gravity water waves to wind excitation. J. Fluid Mech., 35(04):657— 
675, 1969. 

Gottifredi and G. Jameson. The growth of short waves on liquid surfaces under the action of a wind. 
Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences, 319(1538):373-397, 
1970. 

^®W. S. Wilson, M. L. Banner, R. J. Flower, J. A. Michael, and D. G. Wilson. Wind-induced growth of 
mechanically generated water waves. J. Fluid Mech., 58(3):435-460, 1973. 

Mitsuyasu and T. Honda. Wind-induced growth of water waves. J. Fluid Mech., 123:425-442, 1982. 

^^Y.S. Tsai, A.J. Grass, and R.R. Simons. On the spatial linear growth of gravity-capillary water waves 
sheared by a laminar air flow. Physics of Fluids, 17:095101, 2005. 

^^P. Gondret and M. Rabaud. Shear instability of two-fluid parallel flow in a Hele-Shaw cell. Phys. Fluids, 
9:3267-3274, 1997. 

^^P. Gondret, P. Ern, L. Meignin, and M. Rabaud. Experimental evidence of a nonlinear transition from 
convective to absolute instability. Phys. Rev. Lett., 82:1442-1445, 1999. 

^^G. M. Hidy and E. J. Plate. Wind action on water standing in a laboratory channel. J. Fluid Mech., 
26(04) :651-687, 1966. 

^'^G. Caulliez, N. Ricci, and R. Dupont. The generation of the first visible wind waves. Physics of Fluids, 
10(4):757-759, 1998. 

^^G. Caulliez, R. Dupont, and V. Shrira. Turbulence generation in the wind-driven subsurface water flow. 
In Transport at the Air-Sea Interface, pages 103-117. Springer, 2007. 

^^S. Longo, D. Liang, L. Chiapponi, and L. Aguilera Jimenez. Turbulent flow structure in experimental 
laboratory wind-generated gravity waves. Coastal Engineering, 64:1-15, 2012. 

^^G. H. Keulegan. Wind tides in small closed channels. Journal of Research of the National Bureau of 
Standards, 46:358—381, 1951. 

^®J. Francis. Wave motions and the aerodynamic drag on a free oil surface. Phil. Mag., 45:695—702, 1954. 

O Naraigh, P. Spelt, O.K. Matar, and T.A. Zaki. Interfacial instability in turbulent flow over a liquid 
film in a channel. International Journal of Multiphase Flow, 37(7):812-830, 2011. 

“^^J. W. Miles. On the generation of surface waves by shear flows. Part 3. Kelvin-Helmholtz instability. J. 
Fluid Mech., 6(04):583-598, 1959. 

“^^M.-Y. Lin, C.-H. Moeng, W.-T. Tsai, P. P. Sullivan, and S. E. Belcher. Direct numerical simulation of 
wind-wave generation processes. J. Fluid Mech., 616:1-30, 2008. 

“^^F. Moisy, M. Rabaud, and K. Salsac. A synthetic schlieren method for the measurement of the topography 
of a liquid interface. Exp. Fluids, 46:1021-1036, 2009. 

“^^D. Kiefhaber, S. Reith, R. Rocholz, and B. Jahne. High-speed imaging of short wind waves by shape from 
refraction. Journal of the European Optical Society-Rapid publications, 9:14015, 2014. 

‘^'^Sir H. Lamb. Hydrodynamics. Sixth edition, Cambridge University Press, 1995. 

“^^J. C. Padrino and D. D. Joseph. Correction of Lamb’s dissipation calculation for the effects of viscosity 
on capillary-gravity waves. Physics of Fluids (1994-pTesent), 19(8):082105, 2007. 

“^^J. Lighthill. Waves in fluids. Cambridge University Press, Cambridge, 1978. 

“^^H. Schlichtling. Boundary Layer Theory. Springer, 8th edition, 2000. 

“^^S. Longo. Wind-generated water waves in a wind tunnel: Free surface statistics, wind friction and mean 
air flow properties. Coastal Engineering, 61:27—41, 2012. 

“^^A. Zavadsky and L. Shemer. Characterization of turbulent airflow over evolving water-waves in a wind- 
wave tank. Journal of Geophysical Research: Oceans (1978-2012), 117(C11), 2012. 

^^E.J. Plate and G.M. Hidy. Laboratory study of air flowing over a smooth surface onto small water waves. 
Journal of Geophysical Research, 72(18):4627-4641, 1967. 

^^J. Wu. Wind-induced drift currents. J. Fluid Mech., 68(01):49-70, 1975. 

^^Z. Hu, Ch. L Morfey, and N. D. Sandham. Wall pressure and shear stress spectra from direct simulations 
of channel flow. AIAA journal, 44(7):1541—1549, 2006. 

^^J. Jimenez and S. Hoyas. Turbulent fluctuations above the buffer layer of wall-bounded flows. J. Fluid 
Mech., 611:215-236, 2008. 



16 


^'^P. Huerre and M. Rossi. Hydrodynamic instabilities in open flows. COLLECTION ALEA SACLAY 
MONOGRAPHS AND TEXTS IN STATISTICAL PHYSICS, pages 81-294, 1998. 

W. Miles. On the generation of surface waves by shear flows. Part 2. J. Eluid Mech., 6:568-582, 1959. 

J. Plant. A relationship between wind stress and wave slope. Journal of Geophysical Research: Oceans 
(1978-2012), 87(C3):1961-1967, 1982. 

^^T.R. Larson and J.W. Wright. Wind-generated gravity-capillary waves: Laboratory measurements of 
temporal growth rates using microwave backscatter. J. Fluid Mech., 70(03):417—436, 1975. 

^®R. L. Snyder and C. S. Cox. A field study of the wind generation of ocean waves. J. mar. Res., 24:141-178, 
1966. 

Choi and P. Moin. On the space-time characteristics of wall-pressure fluctuations. Physics of Fluids 
A: Fluid Dynamics (1989-1993), 2(8):1450-1460, 1990. 



